New method for measuring longitudinal fluctuations and directed flow in 

ultrarelativistic heavy ion reactions 



L.P. Csernai^^^, G. Eyyubova^ and V.K. Magas'' 
^ Institute of Physics and Technology, University of Bergen, AUegaten 55, 5007 Bergen, Norway 
^ Wigner Research Centre for Physics, 1525 Budapest, Hungary 
^ Department of Physics, University of Oslo, POB 1048 Blindern, 0316 Oslo 
* Departament d'Estructura i Constituents de la Materia, 
Universitat de Barcelona, 08028 Barcelona, Spain 
(Dated: March 3, 2013) 

It has been shown in recent ALICEQLHC measurements that the odd flow harmonics, in par- 
ticular a directed flow, Hi, occurred to be weak and dominated by random fluctuations. In this 
work we propose a new method, which makes the measurements more sensitive to the flow patters 
showing global collective symmetries. We demonstrate how the longitudinal center of mass rapidity 
fluctuations can be identified, and then the collective flow analysis can be performed in the event- 
by-event center of mass frame. Such a method can be very effective in separating the flow patterns 
originating from random fluctuations, and the flow patterns originating from the global symmetry 
of the initial state. 

PACS numbers: 25.75.-q, 25.75.Ld, 25.75.Gz, 25.70.Pq 



Directed flow was the first and most dominant flow pat- 
tern in early relativistic heavy ion studies [1 . It provided 
a solid base for the subsequent more detailed works with 
higher multipole components. With increasing beam en- 
ergy the magnitude of the directed flow decreased as the 
longitudinal momentum increased considerably, and this 
made the measurement of this flow component more dif- 
ficult. 

Recently the strong elliptic flow was demonstrated at 
LHC, exceeding all measurements at lower energies for 
peripheral collisions [5]. At the same time recent mea- 
surements indicate weak directed flow, dominated by 
fluctuations [3 . 

Nevertheless the experience from early studies tells us 
that identification of this flow component is possible and 
it can be separated from a strongly fluctuating back- 
ground. Our goal in this article is to propose a new 
method, which makes the analysis of the flow patterns 
to be more sensitive to global collective symmetries, by 
identifying and then correcting for the longitudinal cen- 
ter of mass (CM) rapidity fluctuations. The proposed 
idea is rather general and can be implemented by all the 
experimental teams, although the particular way of iden- 
tifying the CM rapidity event-by-even will depend on the 
particular set-up. 

Collective flow is parametrized, by the azimuthal angle 
distribution around the beam axis via the expansion 
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where y is the rapidity and pt is the transverse momen- 
tum and (j) is the azimuth angle in the transverse plane 
with respect to impact parameter vector. The functions 
vi{y,pt), V2{y,Pt), ■■■ are the harmonic flow expansion pa- 
rameters. 

Global Collective flow patterns, which follow the global 



symmetries of the reaction Event-by-Event (EbE) pro- 
vide valuable information of the overall pressure and 
transport properties of matter. 

Random fluctuations, especially in the initial state, can 
also lead to flow processes, these however, are not cor- 
related with the global collision geometry and the corre- 
lation of the major axis of the asymmetry arising from 
the fluctuation may have no correlation with the reaction 
plane at all 4j. 

The goal of the present work is to identify and sepa- 
rate the Global Collective dynamical features from ran- 
dom fluctuations. We present some of the steps of such 
reconnaissance for Global Collective flow in high energy 
heavy ion reactions. This is particularly important for 
those harmonic components which are weak and difficult 
to identify. 

The odd harmonics are dominated by fluctuating flow 
components which make it difficult to identify the weaker 
global collective flow components. Just as the transverse 
plane fluctuations, the beam directed fluctuations also 
modify the initial shape |SJ [B] , the tilt of the longitudinal 
axis and similarly the CM position. This issue is not 
discussed in the literature, although it influences all odd 
harmonics in the u„ expansion. 

Recent LHC data indicate that the reaction plane can 
be measured with forward and backward calorimeters 
where the projectile and target spectator residues are de- 
tected providing a reliable detection of the global Event 
reaction Plane EP and the corresponding azimuthal 
angle, '^ep- 

To say something concrete, below we will mostly con- 
centrate on the first flow component, vi, although, as 
it was already mentioned, the method is rather general 
and interesting for all odd components. According to a 
recent estimate the vi flow shows a strong mirror anti- 
symmetric structure as a function of rapidity, but this is 
smoothed out by random fluctuations of the CM motion 
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of the participants. Furthermore it was pointed out that 
fluctuations may cause turbulence i.e. random rotation 
[7]. Although, this was analyzed in the transverse plane, 
similar fluctuating rotation may also appear in the re- 
action plane, which may further soften the directed flow 
peak. 



I. LONGITUDINAL CENTER OF MASS 
FLUCTUATIONS 

The idea is straightforward, if the participant y'~^'^ or 
^CM g^j-Qjigly fluctuating, one can measure it experi- 
mentally, and take the EbE CM into account when the 
odd flow components are evaluated. More than two dozen 
years ago a similar idea was developed, to estimate EbE 
the azimuthal angle of the reaction plane better, and en- 
able the directed flow analysis in low multiplicity and low 
acceptance reactions also [Sj. 

If the acceptance cowers a large fraction of the mo- 
mentum space of emitted particles, the initial CM of the 
system and the final measured CM are nearly identical. 
The final measured CM might deviate from the original 
one if a substantial part of particles are not detected, 
especially if the not-detected particles are not evenly dis- 
tributed in the momentum space. 



A. Participant rapidity from emitted particles 

The total 4-momentum of all measured particles of one 
event is given by 
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where M is the measured multiplicity of the event. P 
can be measured accurately for the pseudo-rapidity ac- 
ceptance range of the detector, I77I < rjmax- The arising 
Center of Mass (CM.) rapidity is 
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If we do not have a good mass resolution the deter- 
mination of Ei, and therefore of E, may become prob- 
lematic, so only the pseudo-rapidity of the CM can be 
determined: 
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If the rapidity acceptance of the detector is limited, 
then the measured longitudinal single particle momenta 
are also constrained. On the other side, the transverse 
momenta are not constrained. For example at ALICE 
TPC detector the rapidity acceptance of the detector is 
limited to \y\ < 0.9, what means that the measured longi- 
tudinal nucleon momentum can not be more than p| | = 1 



GeV/c. At the same time the nucleon transverse mo- 
mentum distribution peaks at pt = 1.2 GeV/c [9], and 
a large part of the distribution extends to a few GeV/c. 
Thus this constraint of the acceptance would result in 
underestimating the CM rapidity. 

The pseudorapidity distribution for peripheral events 
was analyzed in a simple theoretical few source model |10| 
and for a detector with an acceptance of \y\ < 0.9 the re- 
duction in CM pseudorapidity fluctuation was estimated 
to be \t]^'^''\ < 0.2 . This is an estimate relevant for the 
ALICE TPC detector. In other detectors the acceptance 
might be wider, however, due to geometric limitations 
such an underestimation of l??*^*^! may sill occur. Thus, 
we attempt to estimate 77*^*^ in other way. 



B. Participant rapidity from spectators 

Let us consider that we have three subsystems: (A) 
projectile spectators, (B) participants, and (C) target 
spectators. We can measure the energies of A and C, 
Ea and Ec-, in the respective Zero Degree Calorimeters 
(ZDC): ZDCa and ZDCc- Then the energy and momen- 
tum conservation gives 

Eb = Ab mBj_ cosh(?/^) = Etot - Ea- Ec , 
Mb ^ Ab mB± smh{y^) ^ ~{Ma + Mc) (5) 

For example, at the present LHC Pb-|-Pb reaction with 
energy per Nucleon ejv = 1.38 TeV/nucleon, the beam 
rapidity is yo — 7.986 and 

Etot = 2ApbmN cosh(yo) , 

where mjv = 938.8 MeV/c^ 
Furthermore the equations 

Ea = Ap TUN cosh(?/o), 
Ec At tun cosh(-yo), 

give the spectator numbers, Ap and At, and 

Ma = Ap TTiN sinh(yo), 
Mc = Atttin sinh(-yo), 

as well as the mass number of subsystem B: 
Ab = 2Apb - Ap-At . 

Thus for an event, dividing the second of eq. ^ by the 
first we can determine the rapidity of subsystem B, which 
should be close to the rapidity of the participant system. 



(6) 



Our system B includes high energy "pre-equilibrium" 
particles, which are not detected by the ZDCs, and do 
not form a locally equilibrated system. To separate these 
two components from one-another would need more in- 
formation, and a quantitative definition. 
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The part of the neutrons of the colliding nuclei, which 
are originating from the spectators and thus flying with 
beam rapidity (and energy) can be detected by the for- 
ward and backward neutron Zero Degree Calorimeters. 
These ZDCs return the total energy of the measured neu- 
trons. The ZDCs are between the two (Projectile and 
Target) beam pipes. The colliding nuclei propagate be- 
fore and after passing through the point of intersection in 
these beam pipes. Magnetic fields deflect these and other 
charged fragments, so that only the uncharged neutrons 
end up in the ZDCs. 

In central collisions the spectators contain very few, 
mainly single nucleons, and thus ZDC energies tend to 
zero on both sides, A and C. 

At higher impact parameters, the spectators will be- 
come more massive and the energy deposited in the ZDCs 
is increasing. At peripheral collisions, two residue spec- 
tators are expected at opposite sides of the participant 
zone. These may contain single protons and neutrons as 
well as bound nuclear fragments. Of these only the sin- 
gle neutrons reach the ZDCs, because all charged frag- 
ments are deflected away from the the joint beam direc- 
tion (where the two colliding beams interact, and then 
are separated into the two beam pipes). 

In the LHC Pb-|-Pb reactions at 1.38 GeV/nucleon 
beam energy in each beam, results in a total neutron 
energy including both of the colliding Pb nuclei j ~ 
2Npb TTiN cosh(?;o) = 348 TeV, where 2Npb = 2 x 126 is 
the total neutron number in the collision. The neutron 
ZDCs measure only single neutrons. At zero impact pa- 
rameter we do not expect any spectators, while at the 
maximum impact parameter, b = 2Rpb , we do not ex- 
pect any single neutrons, as the two spectators are two 
whole Lead nuclei. At very small impact parameters the 
number of single neutrons are proportional to the num- 
ber of spectators Nn{b) — {N/A)Nspect{b), where the 
nucleon number of the spectators in a collision at im- 
pact parameter 6, can be obtained from the intersection 
geometry (like in the Glauber model). In the rest of 
the impact parameter range the spectators undergo nu- 
clear multi-fragmentation, which is a well studied area 
both theoretically and experimentally [TT] . These multi- 
fragmentation studies enable us to estimate the number 
of single neutrons for a given impact parameter b. 

We can assume that at an intermediate impact param- 
eter where the single neutron number of spectators is near 
to its maximum, we have less than half of the total neu- 
tron number of the spectators. Thus, we can assume the 
maximum of total ZDC energy of about El^f — 160 TeV 
In other words with this estimate from of the 2 x 126 
available neutrons, at most about 2 x 58 single neutrons 
can reach the ZDCs. With increasing impact parameter 
the number of bound nuclear fragments in the spectators 
increases, and thus these events are populating again the 
lower EzDC domains. 

This feature can be checked by plotting E^^^ — 
EzDCa + EzDCc as a function of the total charged mul- 
tiplicity, Nch measured, in the TPC. As discussed above 



the measured E^^^ first increases with decreasing Nch^ 
then reaches a maximum at iV^"* 100 — 200 and then 
it decreases again. This enables us to separate the dis- 
cussion of medium-peripheral collisions with multiplicity 
below N^"* and extreme-peripheral collisions with mul- 
tiplicity higher than N^"*. And thus we can separate 
the correlations between Ez oca £^nd Ezdcc for different 
multiplicity or centrality bins. 

For example this feature of the ZDCs for the ALICE 
ZDCs is discussed in ref. fT^ in an early simulation, 
where incomplete spectator fragmentation is taken into 
account, and HIJING is used as event generator. For the 
estimates a model based on the ALADIN Au-Au data and 
deuterons are estimated from NA49 data. The correla- 
tion between the energy in the Zero Degree Calorimeters, 
versus the number of single spectator neutrons, based on 
theoretical model estimates raised monotonically up to 
about 100 neutrons and then it dropped rapidly. This 
rapid change was clearly attributed to the fact that at 
higher impact parameters the measured neutron number 
does not include neutrons in bound fragments, and thus 
the ZDC energy peaks at a critical neutron multiplicity 

The charged particle multiplicity does not provide a 
good reference measure for peripheral collisions, due to 
small multiplicities and therefore large fluctuations. As it 
is introduced in ref. [2] we can have better visualization 
of the spectator energy, if we plot the total ZDC energy, 
= £'2 + E^, against the event centrality percentage. 

The different presentations of these data are connected 
by the definition of centrality bins, which depend on the 
detector acceptance and the treatment of the lowest and 
highest multiplicities. The connection is presented in Fig. 
1 of ref. !2^. 

Now the above considerations and evaluation of ycM 
should be reconsidered due to the fact that the ZDCs 
measure single neutrons only and the number of single 
neutrons depends in a special way from the impact pa- 
rameter due to the formation of nuclear fragments in the 
spectators or due to the incomplete dissociation of the 
bound spectator nuclear fragment in to single nucleons. 

We can do a quantitative estimate of ycM, if we con- 
clude on the total spectator energy carried by all neutral 
and charged fragments together. Such a correction can 
be based on the measured Nch and N^J^^'^. 

We may attribute the charged particle multiplicities 
or the centrality percentage to the obtained multiplic- 
ity estimated in a fluid dynamical model calculation, 
see Figure [T] which is based on the calculations dis- 
cussed in details in The impact parameters of the 
fluid dynamical model calculation can be matched well 
to the centrality percentages. The correspondence be- 
tween impact parameter and event centrality percent- 
age is impressive for semi-peripheral reactions: b = 
0,2,0.3,0.4,0.5,0.6,0.76™a:r correspond to 5, 13, 21, 29, 
37, 45 % centrality percentage respectively. 

We should correct the preliminary estimate, eqs. (5|6 1, 
for the CM rapidity. We have to correct for the fact that 
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FIG. 1. Simulation of the Pb+Pb collision at LHC with ejv = 1.38 TeV/nucleon [5 . The calculated charged particle 
multiplicity, Nch, as a function of FO time (assuming a tpo = const. FO hyper-surface), for different impact parameters, 
b — 0.0, 0.1, 0.2, ...0.7&mti2!. The indicated (bO, bl, ... b7) FO times for different impact parameters reproduce the measured 
charged particle multiplicities, Nch, in the corresponding centrality bins. The visible fluctuations arise from the feature of the 
PIC method [S], that the volume increases by one cell when a marker particle crosses the boundary. Thus at the initial state 
with relatively few cells and large relative surface, this leads to fluctuations. 



nuclear fragments and the energy carried by these are 
not measured by the ZDCs, and thus our estimates for 
Ea,Ec and thus Eb should be modified and should be 
calculated based on the neutron energies measured in the 
ZDCs. 



1. Neutrons in nuclear fragments 

The simplest approximation is that the total spec- 
tator energies, Ea and Ec are related the same way to 
the measured ZDC energies: 

Ea = iA/N)E^ , Eb = {A/N)El , 

where A/N is the mass number to neutron number ra- 
tio in colliding nuclei (for simplicity we consider sym- 
metric collisions). This approximation is satisfactory for 
central collisions and for small impact parameters (i.e. 
h < 0.3bmax so, relatively large charged multiplicities). 
Due to the presence of bound nuclear fragments in the 
spectators, we need a better approximation at inter- 
mediate or higher impact parameters. 

Based on the initial state model [13], used in the 3 di- 
mensional fluid dynamical calculations [SI [T3| from the 
geometrical overlap of the colliding nuclei, we get the 



number of participant nuclcons for each impact param- 
eter. On the projectile (or target) side 126/208 part of 
these are participant neutrons, Pn(b), shown by dashed 
blue line in Figure [2j The remaining neutrons are in spec- 
tators, shown by the dotted green line. As in ref. [T3j, we 
assume that the N/Z ratio is homogeneous in the whole 
initial state system. 

Not all spectator neutrons are single neutrons; some 
of these are in composite nuclear fragments. The total 
energy of single neutron spectators are measured in the 
ZDCs, and thus their number can be easily obtained as 
the beam energy per nucleon is known. So, we can get the 
number of single neutrons, Nn (6) , as a function of the im- 
pact parameter 6 in a straightforward way from the total 
measured ZDC neutron energy EzDCa + Ezdcc — E^^f 
as a function of the impact parameter or the central- 
ity percentage. The asymptotic behaviour of this de- 
pendence is not straightforward, i?"^^ — > both if 
b/bmax — ^ or 1. Based on nuclear multi-fragmentation 
studies mentioned above 11 , we can estimate the num- 
ber of single spectator nucleons, iV„(6), as given in Table 
m This estimate indicates that the number of detected 
single neutrons should have a maximum of about 30% at 
^ 50% centrality percentage. 

The difference between the average and maximum for 
the single neutron numbers leads to an estimate of the 
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TABLE I. The number of single spectator neutrons as func- 
tion of centrality bins and the corresponding impact parame- 
ters as estimated based on nuclear multi-fragmentation model 
|11| . The initial geometrical spectator numbers corresponding 
to a given impact parameter are also given. 

systematic error of the method, for estimating the total 
spectator energy and momentum from the single neu- 
trons. 

Participant and Spectator Neutrons 
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FIG. 2. (color online) The number of participant neutrons 
from the projectile or from the target (dashed blue line) and 
the corresponding number of spectator neutrons (dotted green 
line) in one of the spectators (forward or backward) obtained 
in the initial state calculation [13] for Pb-I-Pb collision.. The 
neutron distribution is assumed to be homogeneous in the 
system, N/A — 126/208. At large impact parameters part 
of spectator neutrons remain in nuclear fragments, which are 
charged and so do not reach the neutron ZDCs. Based on the 
FD estimates we relate impact parameter with the centrality 
percentage, and based on nuclear multi-fragmentation studies 
|11| . see Table [l] we estimate the number of single nucleons, 
which reach the neutron ZDCs (magenta dashed double dot- 
ted line). 

Part of these spectator neutrons are inside compos- 
ite nuclear fragments, especially at large impact param- 
eters. These are charged, and due to the magnetic fields 
directing the beam, these do not reach the central neu- 
tron ZDCs nor the proton ZDCs. As discussed above 
his reduces the energy detected in the ZDCs to about 
130-140 TeV, and the maximum numbers of detected sin- 



gle neutrons to about 47-51. This maximum appears at 
the centrality percentage of ~ 50%, which corresponds 
to & = 0.75bmax according to our fluid dynamical model 
estimates. 

Thus we estimate number of single neutrons in the 
spectators for each impact parameter. These estimates 
have some systematic error due the unavoidable theoret- 
ical input, and at very central collisions as well as for 
extreme peripheral collisions (where there are no data 
and the multiplicity is so low that the fluid dynamical 
approximation is not applicable). The multifragmenta- 
tion process may have a fluctuation itself, which may 
interfere with the Ea Ec asymmetry, and lead to some 
systematic error. 

In a given experimental setup one can measure directly 
7V„(6), via measuring the average of the total energy in 
the two ZDCs as function of the impact parameter, 6, 
and dividing it by the beam energy per nucleon. This 
direct measurement of iV„(&), should be preferred, as it 
eliminates large part of the uncertainties of the theoret- 
ical approach arising from directly emitted high energy 
neutrons, and special emission evaporation mechanisms 
at the spectator/participant boundary with large shear. 

Now we have to estimate, at a given centrality per- 
centage (or impact parameter, 6), the total energy of the 
spectator residues (including protons and charged frag- 
ments), Ea and Ec, from the energy of measured single 
neutrons in the ZDCs, E^ and Sp". From the measured 
spectator neutron energies, E^ and E(V-, we get the cor- 
responding single neutron numbers, = E'^^/epf and 
N^'^ = E^^/eN- Multiplying these numbers by the cor- 
rection factor, Pn{b)/Nn{b), we get the spectator num- 
bers, Ap^T = N^(jPn{b)/Nn{h), and spectator energies, 
including the contributions of single protons and of all 
nucleons bound in composite nuclear fragments: 

EA{b) = {A/N)ETPn{b)/N^{b) 

Ec{b) = (A/iV)£;j"P„(6)/7V„(6) . (7) 

This yields the corresponding spectator momenta, 
M^,M(7, and we can get the event by event CM. ra- 
pidity as in eq. ^ 

y^^b) ^y^^ artanh J^-^ + Me \ yCM^^^^^ 

(8) 

where the last term is added to correct for eventual de- 
tector asymmetry, which is measurable, for all events of 
the sample for a given multiplicity percentage bin. Now, 
EatEcAp^Ax are estimated and in the estimate we 
used the average Nn{b), based on the estimated or even- 
tually measured ZDC energies. 

The correction increases for increasing impact param- 
eter, which leads to increased estimates for ycM fluctu- 
ations, on the other hand the correction is also increases 
and this leads to the possibilities of larger systematic er- 
rors. 

At extreme peripheral reactions spectators appear in 
the form of bound, and so charged, nuclear fragments. 



6 



which are diverted by the magnetic field and are not 
detected. These provide more room for large CM. ra- 
pidity fluctuation, but unfortunately neither the ALICE 
TPC nor the ZDCs can measure these reactions with the 
needed acceptance. 

The above suggested analysis, by using the ZDC and 
charged multiplicity information together may still pro- 
vide the best estimate for the longitudinal fluctuations. 



C. Participant rapidity fluctuations 

In central collisions the ellipticity from random fluc- 
tuations can easily supersede the one from the global 
symmetry. This is indicated by the recent flow measure- 
ments, where the axes of different higher harmonics are 
uncorrelated, indicating that the source is random fluctu- 
ations, which are not correlated with the global reaction 
plane. 

The global asymmetry may become dominant at higher 
impact parameters, i.e. b > 0.6 — 0.8b„iax- For example, 
using the observed multiplicity in the ALICE TPC the 
centrality bins were defined as shown in Ref. These 
centrality bins can be assigned to impact parameter in- 
tervals, based on theoretical and geometrical model as- 
sumptions. 

For example, in Ref. [S] simulations of the Pb+Pb 
collisions at LHC have been performed. The compari- 
son of the multiplicity estimates of this model provide a 
correspondence between the impact parameter and the 
centrality percentage, see Fig. [l] 

For the azimuthal flow asymmetries arising from ran- 
dom fluctuations, we cannot expect a correlation be- 
tween the reaction plane (as measured by the spectator 
residues), and the asymmetry observed in the TPC. This 
also applies to the center of mass of the observed par- 
ticles, because the longitudinal fluctuations are also not 
correlated with the spectator residues. Thus, we can ex- 
pect a correlation mainly at peripheral collisions, where 
the ellipticity arising from global symmetry may be sig- 
nificant and may exceed the one caused by random fluc- 
tuations. Based on the above discussed correspondence 
between the multiplicity and the impact parameter, we 
can expect to see global effects at 50% centrality bin or 
higher, i.e. in rather peripheral collisions. 

The collective global flow components appear at mod- 
erate momenta, thus the CM determination should be 
sufficient for this component of the flow. In the previ- 
ous subsection two methods were presented to estimate 
the CM. rapidity. Using the TPC with the restricted 
in pseudorapidity acceptance limits the observed longi- 
tudinal momenta, while the transverse momenta are not 
constrained. Due to this the TPC is expected to under- 
estimate the CM. rapidity fluctuations. 

The other method based on the single neutrons in the 
ZDCs is not constrained the same way. On the other 
hand one has to estimate the total spectator energy and 
momentum from the single neutrons, and this step re- 



quires some theoretical estimate which introduces a sys- 
tematic error. This is illustrated by the fact that the 
correction factor, P„(fe)/A^„(6) can be based on the av- 
erage single neutron multiplicity. Nn{b) can also be ob- 
tained from the theoretical estimates or the measured 
data, where the total single neutron energy, measured in 
the two ZDCs, is divided by the beam energy per single 
neutrons, ejy. 

The EbE data would show large fluctuations within 
each event among the emitted particles, exceeding the 
CM. fluctuations considerably. We expect that com- 
pared to these fluctuations the CM. fluctuations are 
small. 



D. Correlation between the TPC and ZDC CM. 
rapidities 

Before we start to study the correlations we have to 
remove the average CM. rapidity shift from the data 
as this arises from asymmetries in detector acceptance. 
This is probably negligible for the TPC but expected to 
be significant for the ZDC data. 

The azimuthal correlations of global collective flow are 
in correlation with the reaction plane, also the longitudi- 
nal CM. rapidity should be globally defined for the flow 
pattern. In such a situation one would expect that the 
TPC and ZDC rapidities should strongly correlate. 

On the other hand if the longitudinal and azimuthal 
correlations arise from random fluctuations, which are 
not correlated with the global symmetry axes then such 
a correlation for CM. rapidities measured at different 
rapidity ranges should not correlate. 

We anticipate that at central and semi-central colli- 
sions there is no significant correlation or anti-correlation 
among the CM. rapidities shown by the particles de- 
tected in the TPC or the ZDCs. In these collisions the 
correlations are arising dominantly from random fluctu- 
ations. 

For peripheral collisions of 50 % centrality percent- 
age or more a weak positive correlation is anticipated. 
Here the global collective asymmetry (eccentricity) is so 
strong, that the global correlations start to become com- 
petitive. 

II. ADJUSTMENT OF THE CENTER OF MASS 

Substantial initial state rapidity fluctuations will av- 
erage out all flow structures around the CM when the 
measurements are assuming that the CM is identical with 
the pre-collision CM of the given experiment (i.e. fixed 
to the Laboratory frame in a colliding beam experiment 
with a symmetric, A-l-A, collision). Odd components of 
global collective flow patterns, {vi, vs, ...) are Mirror 
Asymmetric (MA) around the real participant CM, so 
these are severely effected by the EbE CM fluctuations 

mm- 
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Based on the above results, we suggest to use the CM 
rapidity, EbE, determined from the ZDC data. Let us 
assume that the CM rapidity is measured for each Event, 

where P = "Yli^^iPv the total measured 4-momentum 
of event E. Here we consider that all corrections are 
taken into account as described in eq. (|8]). Let us then 
shift each event to its own CM by the measured y^*^ so 
that each particle rapidity will be moved to 



(9) 



This transformation will not effect the azimuth angle of 
the emitted particles, nor m±^i, however, the longi- 
tudinal momentum and the energy will change to, 



El = m±,i cosh(yi 



Consequently the EbE Vn{p±) will not change after ad- 
justing the center of mass rapidity (if we integrate over 
the full rapidity range). Instead in a given event, E, this 
adjustment will contribute to a re-distribution to 
of at another set of rapidities, y-. Consequently the de- 
tector acceptance boundaries will also change EbE, and 
so one will have a continuous rapidity coverage even if the 
detectors have some rapidity acceptance gaps. Then, the 
flow harmonics can be determined EbE averaging over 
all measured particles in the Event 



Vn{y',P±)E = ( C0S[71( 



EP)\ IE 



(10) 



and then one can make an average over all events in a 
centrality bin: 



Vn{y\Pl_) = ( Vn{y',P±)E ) 



(11) 



If the centrality bin is wide this two step averaging is 
essential [15^. 

As usual, one can also correct for fluctuations of the ob- 
served Event Plane azimuth angle, compared to the 
real (pre-coUision, but not measurable) Reaction Plane, 
\E'_R, as 7^ "^EP, and so the observed must be 
corrected by dividing by the resolution of the event plane 

m- 

So, using eqs. (8|9) we can perform this rapidity shift 
event-by-event. To eliminate a relatively large rapidity 
binning of the sample the events could be distributed 
uniformly and randomly inside each rapidity bin. Then 
the pseudo-rapidity distribution of Vi , as measured in the 
EbE CM frame is expected to show a clear ±y asymmetry 
as the longitudinal, y^^^ fluctuations are removed from 
the data. The distribution this way may show a distinctly 
antisymmetric distribution in rapidity, as required by the 
global symmetry of the event. 



In case of ALICE, for example, the TPC data are avail- 
able in terms of pseudorapidity and integrated over the 
p± acceptance of the TPC. In such a situation instead of 
eq. ([9]), we can use the approximation that the pseudo- 
rapidity, EbE, is shifted by the CM rapidity: 



vl = n.-yr\b) 



(12) 



So, this way the collective part of the flow may be 
identified. 

We can quantify the identification of the collective 
symmetric versus the random fluctuating contribution to 
the flow the following way. We can evaluate the odd and 
even components of the flow |3]: 



vn{y',p±) 
vT'^iy'.P^) 



K(y',?3_L) + vn{-y' ,piy\ 
K(y',?3_L) - Vn{-y',pi^)] ■ 



We can use here also the approximation eq. ( 12 ) to eval- 
uate the symmetries of the data. 

Initial states with global collective symmetry are fully 
Mirror Asymmetric (MA) states, while fluctuating states 
may have both MA and Mirror Symmetric (MS) com- 
ponents. Then the two, odd and even components for 
vi{ri') should show a distinct difference. The original 
vi{ri') component is expected to be very close to the odd 
component, while the even component should be much 
smaller in the central rapidity range. This would indi- 
cate that vi{ri') is MA to a large extent, and the MS 
component is consistent with zero. 

At zero impact parameter global symmetry does not 
result in any global azimuthal asymmetry, so all har- 
monic coefficients must vanish from global collective ori- 
gin and only random initial state fluctuations as well 
as dynamical random fluctuations developing during the 
collision may lead to azimuthal fluctuations. The recent 
ALICE analysis of azimuthal flow asymmetries central 
collisions [3] , could identify flow harmonics from vi to 
and shows a maximum for 



III. CONCLUSIONS 

We presented a method how to analyze the directed 
flow data by considering the EbE longitudinal rapidity 
fluctuation of the participants. The participant CM ra- 
pidity can be estimated both from the TPC and the ZDC 
data. The TPC data are constrained to the pseudora- 
pidity acceptance range, which impairs the CM rapidity 
estimate for peripheral reactions. To obtain an estimate 
for peripheral reactions from the ZDC data we estimate 
the forward and backward spectator energies from the 
measured single neutron energies. For peripheral reac- 
tions, where the azimuthal asymmetry of overall global 
collective origin is expected to exceed those, which origi- 
nate from random initial state fluctuations, we expect to 
flnd significant correlation between the TPC and ZDC 
estimates of the CM rapidity of the participants, in the 
central rapidity range. 
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Using the ZDC data, which arc not constrained by the 
limited acceptance, we describe each event in its own CM 
frame, and propose to evaluate the directed flow from 
these shifted data. 



The method of shifting the system origin event by 
event to the rapidity of the participant CM rapidity is 
effective in separating the flow patterns originating from 
random fluctuations, and the flow patterns originating 
from the global symmetry /asymmetry of the initial state. 
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